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A GENERALIZED ORTHOGONAL COORDINATE SYSTEM 
FOR DESCRIBING FAMILIES OF AXISYMMETRIC 
AND TWO-DIMENSIONAL BODIES 

Peter A. Gnoffo 
Langley Research Center 

SUMMARY 

A generalized curvilinear orthogonal coordinate system is presented which can be 
used for approximating various axisymmetric and two-dimensional body shapes of interest 
to aerodynamicists. Such body shapes include spheres, ellipses, spherically capped cones, 
flat -faced cylinders with rounded corners, circular disks, and planetary probe vehicles. 

In addition, a set of transformation equations is developed whereby a uniform velocity 
field approaching a body at any angle of attack can be resolved in the transformed coordi- 
nate system. The expressions for the metric coefficients in this coordinate system are 
derived, and it is found that two of the metric coefficients are identically, equal. The 
Navier-Stokes equations are written in terms of a generalized orthogonal coordinate sys- 
tem to show the resultant complexity of the governing equations. Also outlined is a solu- 
tion procedure to obtain the flow field surrounding a planetary probe vehicle traveling 
at 0® angle of attack at supersonic velocity through a perfect gas. 

INTRODUCTION 

The introduction of vector -processing computer systems has enhanced the feasibility 
of solving the full Navier-Stokes equations for compressible flow over various body shapes. 
However, the efficiency of vector processing is a function of how well (or how simply) a 
problem can be vectorized (ref. 1). In treating the flow over an axisymmetric body, it is 
advantageous if the body and the surrounding flow field can be transformed into a simple 
rectangular coordinate system. The vectorization is further simplified if the boundary 
conditions are transformed to one boundary of a rectangle, and the free-stream conditions 
are transformed to the opposite boundary of the rectangle. 

In this paper an orthogonal coordinate system with these properties is presented for 
analytic body shapes which approximate spherically capped cones, flat-faced cylinders with 
rounded corners, circular disks, and planetary probe vehicles. The orthogonal coordinate 
system will also handle exactly the transformation of bodies with circular or elliptical 
cross sections. 



SYMBOLS 


constants of transformation 
metric coefficient 
total enthalpy 
length of body 

characteristic cone length defined in equation (14) 

slope 

integer 

pressure 

coordinate in transformed space 
radius of curvature 
arc length 
time 

temperature 

velocity along line of constant r 
velocity along line of constant 6 
total velocity 


coordinates in real space 



Y „ maximum body radius 

max 

a angle of attack 

y ratio of specific heats 

9 coordinate in transformed space 

6 „ value of 9 for maximum body radius 

max 

X second coefficient of viscosity 

jLt viscosity 

p density 

T^j shear stress along surface perpendicular to i-axis in direction parallel 

to j-axis 

(f) coordinate in transformed space for axisymmetric case 

ip reference angle for line of constant r in real space 

ojjjWgjCUg velocities in generalized coordinates 
Subscripts: 
b body 

B base 

c cone conditions 

n nose 
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r pertaining to conditions along line of constant 6 and 0 

s shock 

9 pertaining to conditions along line of constant r and <p 

0 pertaining to conditions along line of constant r and 0 

°o free -stream conditions 

DISCUSSION OF THE PROBLEM 

The solution of the Navier-Stokes equations over a two-dimensional or axisym- 
metric body presents many difficulties in the computational formulation of the problem. 

One of the major factors which affects the complexity of the computational formulation is 
the choice of a suitable coordinate system to be used as a basis for differencing the gov- 
erning equations. 

There are basically two options available in the selection of such a coordinate sys- 
tem. The first option is simply to use either a rectangular coordinate system for a two- 
dimensional body or a cylindrical coordinate system for an axisymmetric body. The 
Navier-Stokes equations are in their simplest form (metric coefficients all equal to one) 
in a rectangular coordinate system so there is no need to calculate or store metric coef- 
ficients in a computer program. The Navier-Stokes equations written in a cylindrical 
coordinate system are only slightly more complicated and include a metric coefficient 
which varies as the distance from the axis of symmetry (ref. 2). However, there is no 
guarantee in either system that the body of the problem will pass through node points of 
the finite difference mesh used in the coordinate system (see fig. 1). Mesh points can be 
chosen to lie on the body surface, but the problem becomes more complex because a large 
amount of work is necessary in order to compute mesh points and boundary conditions for 
each different body shape which needs to be examined. This specification of mesh points 
will also force one to use differencing techniques based on unequal spacing at points close 
to the body. On the other hand, if mesh points are not made to lie on the body, some type 
of interpolation scheme has to be written into the program in order to compute boundary 
conditions and property derivatives on the body surface. Another problem with the rec- 
tangular and cylindrical coordinate systems is that there is no simple way of. concentrat- 
ing mesh points near the surface of an arbitrary body. Such a concentration is desirable 
because the viscous effects in a boundary layer lead to large gradients which will be inade- 
quately defined by a coarse mesh. 
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In reference 3, Moretti describes a polar coordinate system (fig. 2) which is more 
realistic for describing two-dimensional and axisymmetric body shapes. This transformed 
coordinate system is generally nonorthogonal but does have the advantage of being able to 
map most arbitrary bodies onto a rectangle. Also mesh points can be concentrated near 
the body. 

The second option available in the selection of a coordinate system is to choose a 
system based on the body itself. Within this option, two subgroups can be defined: an 
orthogonal body-oriented coordinate system and a natural coordinate system. Often an 
orthogonal body-oriented coordinate system can be used to describe bodies composed of 
combinations of spheres, cones, and other analytic body shapes. In this coordinate sys- 
tem, the coordinates of a point are determined by s (the distance along a body surface 
measured from the axis of symmetry) and n (the distance along a line which is normal 
to the body at s). A body-oriented coordinate system can be used, for example, to 
describe a spherically capped cone with a hemispherical afterbody joined smoothly at the 
corners (fig. 3(a)). The major problem with this system in this case is that four different 
sets of metric coefficients would have to be calculated, one set for each analytic section 
of the body. The body curvature, and consequently, one of the metric coefficients and the 
derivative with respect to s of another metric coefficient would not be continuous func- 
tions of s. The matching of flow properties at the boundaries of these various analytic 
sections and the formulation of the differencing procedure at these boundaries would add 
to the complexity of the problem. Also, in a body-oriented coordinate system, any dis- 
continuity in slope on the body makes it impossible to describe a large section of the flow 
field. This situation would occur, for example, in attempting to describe the base flow 
for a sphere-cone configuration (fig. 3(b)). The body-oriented coordinate system would 
also break down when a concavity occurs in the body since the concavity leads to multiple 
specifications of some area in front of the concavity (see fig. 3(c)). 

In a natural coordinate system, the body surface itself forms a boundary of a trans- 
formed computational space. Some obvious examples of this approach are a spherical 
coordinate system to describe flow over a sphere, a cylindrical coordinate system to 
describe flow over a cylinder, and a parabolic coordinate system to describe flow over a 
body with a parabolic cross section (fig. 4 and ref. 4). Metric coefficients for all of these 
coordinate systems are simple to calculate. Moreover, there is little difficulty in con- 
centrating mesh points near the body because the computational mesh system is composed 
of lines parallel to the body which can be concentrated as close to the body as desired. 

One very interesting natural coordinate system which is generally nonorthogonal is 
the "rock transformation" (ref. 5). Transformation functions are generated such that all 
boundaries of a problem are coincident with coordinate lines. The transformation of the 
governing equations results in a complex system of partial differential equations with 
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simple boundary conditions. Other interesting coordinate systems and coordinate trans- 
formations are outlined by Roache in reference 6. 

In this paper, a more general natural orthogonal coordinate system is presented 
which can be used to approximate many different axisymmetric and two-dimensional body 
shapes. An example of this new coordinate system, shown in figure 5, transforms the 
area surroxmding the body in real space into a rectangular computational grid. In this 
grid a finite length of one coordinate axis is the transformed body. As will be shown, the 
coordinate system is able to uniquely describe the area surrounding bodies with concavi- 
ties. As might be expected, the expressions for the metric coefficients become more 
complicated; however, it turns out that two of the metrics are equal to each other. It can 
also be noted that for some of the bodies which can be described by this method a body- 
oriented coordinate system can be used to describe the area surrounding the analytic body. 
However, it is the development of the natural coordinate system that will be discussed in 
the following sections because that system can be used to describe the region surrounding 
all of the bodies presented in this paper. 

GENERAL TRANSFORMATION 


The transformation from the (0,r,</>) domain to the (x,y,z) domain for an axisym- 
metric coordinate system is written as 


N 

x{6,r,(p) = (-B sinh r + C cosh r) cos ® ^ cos(n0) 

n=2 


yiO,r,(p) = 


N 

(B cosh r - C sinh r) sin ^ ^ sin(n0) 

n=2 


cos (p 


z{e,r,(p) = 


N 

(B cosh r - C sinh r) sin 0 + ^ Aj^ e^’^ sin(n0) 

n=2 


sin <p 


( 1 ) 


where N is a positive integer >2 and Aj^, B, and C are arbitrary constants. A two- 
dimensional transformation to the xy-plane is obtained by setting <p = 0. 

Lines of constant r are perpendicular to lines of constant 0 in the xy-plane. 
This relationship can be proved by showing that the slope m^ of lines of constant r is. 
equal to the negative inverse of the slope m^. of lines of constant 0 where such lines 
intersect; that is, 
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9y 

“ 9x 


9y \ (B cosh r - C sinh r) cos 0 + / nA cos(n0) 

\ h 

/r £2S\ N 

It (b sinh r - C cosh r) sin ^ ^ sin(n0) 


(2) 


m_ = 




N 

(B sinh r - C cosh r) sin 0 + ^ sin(n0) 

n=2 




dr 


N 

1 9 (-B cosh r + C sinh r) cos 0 ■ ^ cos(n0) 

n=2 


(3) 


By inspection it is seen that m^. 



Lines of constant r are transformed to circles in the xy-plane in the limit as 
r approaches negative infinity. This transformation can be verified by observing cer- 
tain properties of the transformation. Terms involving the factor e^^^, will vanish 
as r approaches negative infinity. Also note that -sinh r approaches cosh r as 
r approaches negative infinity. Therefore, for large negative r, equations (1) can be 
written approximately as 


x(0,r) a Qb + C) cosh 0 cos 0 




y(0,r) s Qb + C) cosh 0 sin 0 


(r « 0) (4) 


Equations (4), which are valid approximations to the transformation for large nega- 
tive r, are the equation of a circle in the xy-plane with a radius equal to (B + C) cosh r. 
Also note that lines of constant 0 become straight lines with slope mj. equal to tan 0 
in the region where equations (4) are valid. 

The line segment r = 0, 0 ^ 9 < 2n is transformed into a two-dimensional body 
in the xy-plane and becomes one boundary of the computational space in the 0r-plane. An 
axisymmetric body is defined by mapping the line segment r = 0, 0 ^ 9 = tj to the 
xy-plane and rotating the image around the x-axis 2t! rad. Thus, 
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Two dimensions 


N 

= C cos 0 - ^ A„ cos(n0) 
n=2 

N 

= B sin 0 + ^ Aj^ sin(n0) 

n=2 j 


(0 £ 0 < 2n) (5) 


Axisymmetric 

N 

X|^ = C cos 0 - ^ An cos(n0) 
n=2 


yb = 


B 


N 

sin 0 + An sin(n0) 


n=2 


cos 




^b = 


N 

B sin S + ^ sin(n0) 
n=2 


sin 0 


/O i 0 ^ 7T \ 

\0 £ 0 < 277 / 


( 6 ) 


The lines 0 = 0 and 6 = 2tt will therefore form two more boundaries of the com- 
putational space for a two-dimensional problem while the 0 = 0-plane, 0 = jr-plane and 
the 0 = 0-plane, 0 = 27 r-plane are the boundaries in computational space for an axi- 
symmetric body. The final boundary of computational space is specified as r = r^. The 
magnitude of r^ must be chosen large enough so that approximate boundary conditions 
can be applied. 

It should be noted also that this transformation can be written as a conformal map- 
ping for two dimensions. By setting w = x(0,r) + iy(0,r), E = (B + C)/2, and 
Aj^ = (B - C)/2 and by simplifying the result, one obtains 


N 

w = Ee^^ - ^ 


-inz 


n=l 
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where 


z = 0 + ir 


and 



The problem which now remains is to pick the constants B, and C to 

describe the body of interest. 


DETERMINATION OF BODY SHAPE 


At first it appears that it might be possible to determine from equations (5) the con- 
stants An, B, and C for any body shape by using a Fourier transformation. However, 
once the constants for y{9) have been specified, all but one of the constants for x( 0 ) 
are also specified. Therefore, it should not be expected that any arbitrary body shape can 
be produced by some systematic approach to determining the constants in equations ( 5 ). 
Also, once the body shape and the computational space are completely defined, the metric 
coefficients will have to be calculated. The calculation of these metrics (to be described 
later) becomes extremely tedious if N is large. Therefore, the approach taken in 
this section is to see the kind of body shapes that result by defining a small number of 

constants. ,j 

'1 

Families of Ellipses 

For all An = 0 the transformation from the 0r-plane to the xy-plane results in 
families of ellipses. For B = C the transformation results in a circle of radius B 
(fig. 6 (a)). For B 7 ^ C the transformation results in an ellipse (fig. 6 (b)) of eccen- 
tricity e, where e is defined by the following equation (ref. 7): 


e = 



( 7 ) 


Families of Triangular Shapes 

Families of triangular shapes occur when A 2 7 ^ 0 and A ^^>2 = 7'he radius of 
curvature on the body can be written as follows (ref. 7): 
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,3/2 


3/2 


R = + 



ydx 




ax 3^y _ ay a^x 
90 qq2 dO 3q2 


( 8 ) 


The radius of curvature becomes infinite when the denominator of the last expres- 
sion in equation (8) equals zero and the numerator is finite. For B = C = 1 and 
Ag = 1/4, one can substitute into equation (8) to find that an infinite radius of curvature 
occurs for 0=0, 27 t/ 3, and 4^3 (fig. 7(a)). For A 2 > 1/4 or A 2 < 1/4 the body- 
shape change appears as is shown in figure 7. Care must be taken in that if A 2 is too 
large the body contour will intersect itself (fig. 7(c)). This condition is unacceptable for 
defining a real body shape in the xy -plane. 

Also, in this three parameter system where A 2 , B, and C are free parameters, 
there is the option of specifying three independent geometric conditions for a particular 
body type, and then it is possible to solve for A 2 , B, and C. This approach was taken 
to define the bodies in figure 8 by specifying a nose radius of curvature at 0 = n, 
an infinite radius of curvature at 0 = 0, and a maximum body radius 

Families of Rectangular Shapes 

Families of rectangular shapes occur when A 2 = 0, Ag / 0, and A^^^g = 0. If 
equations (5) are substituted into equation (8), then for B = C = 1 an infinite radius of 
curvature exists at 0 = 0, ir/2, ir, and 3 t 7/2 for Ag = 1/9. The resulting body shape 
(fig. 9) approximates the cross section of a flat -faced cylinder with rounded corners. As 
with the triangular shapes, three independent geometric conditions can be specified in 
order to solve for Ag, B, and C. By specifying an infinite radius of curvature at 
0 = 0, a total body length and a maximum body radius Y„„„, the body shapes in fig- 
ure 10 were constructed. Note that the body outlined in figure 10(b) can be used to 
approximate the cross section of a flat, circular disk. 


Families of (N+1) -Sided Shapes 

For B = C = 1 and Aj^ = with all other Aj^ equal to zero, a family of 


(N + l)-sided shapes is produced. This can be shown as follows. Let 


1 ^ 

X. (0) = cos 0 cos(N0) 

^ xt2 


n2 

y. (0) = sin 0 + ^ sin(N0) 




( 9 ) 
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Substituting the above expression into equation (8) yields 


R = 


{■ 


-sin 6 + ^ sin(N0) 


cos 0 + ^ cos(N0) 


v3/2 


-sin 0 ^ sin(N0) 

N 


-sin 0 - sin(N0) 


cos 0 -I- ^ cos(N0) 
N 


-cos 0 + cos(N0) 


( 10 ) 


The denominator of equation (10) represented by D reduces to 


D = (1 - - 1 (1 - cos 
N 


gN f l)^j 


If Q = denominator becomes 

N + 1’ 


D = |l ■ " cos(2k77)J = 0 


( 11 ) 


( 12 ) 


for k equal to 0 or an integer. Thus for N = 3 points of zero curvature (infinite 

2k 71 

radius of curvature) occur at 0 = for k = 0, 1, 2, and 3. For N = 4 points of zero 

4 

2k?7 2k?7 

curvature are produced when 0 = , and for N = 5 they are produced when 0 = 

5 o 

Figures 11 and 12 show the resulting transformation into the xy-plane and demonstrate 
these conditions for N = 4 and N = 5. 


Approximation to Spherically Capped Cones 

A good approximation to spherically capped conical bodies can be obtained by 
including one extra parameter Ag into the family of triangular shapes described in a 
previous section. Thus, 


x^j(0) = C cos 0 - Ag cos(30) - Ag cos(20) 


yjjj(0) = B sin 0 + Ag sin(30) + Ag sin(20) 


(13) 


The constants Ag, Ag, B, and C are determined by simultaneously solving the 
equations specifying basic geometric parameters of a spherically capped cone, namely, 
0g, Yj^^,.and Rj^ (fig. 13). The base of the sphere cone at 0 = 0 is specified as 
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having zero curvature. The radius of curvature of the nose at 6 = n 
The value of for the maximum body radius is determined by setting 

The distance along the symmetry line defined by the equation 


is specified. 
de(^max) = °- 


^c =^b(®max) - 

is determined from the sphere-cone parameters by 


(14) 


^c = 


^max ^c-^n(^-"^"^c) 
sin 


(15) 


The maximum radius of the body is set equal to the base radius of the equivalent 

sphere cone. These equations can now be written as 


C - 9Ag - 4A2 = 0 (16) 

— B — ^Ag + 2A2 
C — QAg + ^^2 



= R. 


(17) 


3A3Cos(39^ax)+2A2Cos(2e„^)tBcos9^^ = 0 (18) 

-A3Cos(3emax) -A2“®(2«max) +'^“®%ax-A3+A2 + '= = 'c 

*3 =‘‘'(39max) + *2 ^‘"(^Vax) i- ® “max ' ^max 

. Equations (16) to (20) can be solved numerically using Newton's method to obtain A 2 , 
Ag, B, C, and Specifications of various spherically capped cones and the equiva- 

lent transformed bodies appear in figures 14 to 18. True sphere cones are plotted over 
these bodies to demonstrate the agreement of the analytic approximation to the true body 
shape. For a ratio of nose radius to base radius equal to 1/2, fairly good agreement is 
obtained in the cone-angle range from 55° to 60°. The approximation becomes less accu- 
rate as the cone angle decreases. The forebodies of the analytic bodies for these low cone 
angles approximate hyperbolas better than sphere cones. A hyperbola of equivalent nose 
radius and asymptotic angle is plotted in figure 18. It is again cautioned that completely 
arbitrary specification must be avoided. For example, specifications for yery wide angle 
bodies with small ratio of nose radius to base radius wilT result in bodies which loop 
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around and intersect themselves (fig, 19). Once a specific sphere-cone configuration is 
picked and the parameters Ag, Ag, B, and C are obtained, the transformed body 
should be plotted to determine whether or not agreement to the desired shape is 
satisfactory. 

It should be noted that lines of constant r are concentrated near the body surface 
and lines of constant 0 are concentrated in regions of sharp curvature. (See, for 
example, fig. 20.) This natural concentration of nodal points near the body surface should 
prove advantageous when solving the flow -field equations around the body using differ- 
encing techniques because the largest gradients in the flow-field variables typically occur 
near the body and at the corners of the body where severe expansion phenomena occur. 

-^Planetary Probe Vehicles 

A commonly considered configuration for planetary probe vehicles is a spherically 
capped conical forebody with a hemispherical afterbody. The forebody and the afterbody 
are blended together through a rounded curve at the shoulder as, in figure 21. This con- 
figimation can be approximated by including one additional term A^ to the sphere -cone 
approximation terms. Thus, 


Xb(0) = -A^ cos(40) - Ag cos(30) - Ag cos(20) + C cos 9 


= ^4 sin(49) + Ag sin(39) + Ag sin(20) + B sin 9 J 


( 21 ) 


As with the sphere- cone approximation, the constants in equations (21) can be deter- 
mined by specifying the basic geometric parameters of a planetary probe vehicle. It is 
possible to solve for A2, A3, A4, B, and C by defining a nose radius of curvature 
Rjj, a base radius of curvature Rg, a vehicle length a maximum body radius Y^ax* 
and the length from the nose to the location of the maximum body radius on the sym- 
metry line l(>. To specify Y^ax H Is also necessary to determine 0max such that 

y(®max) = ^max ^(®max) = These specifications result in the following 

equations: 


-(B + 4 A 4 + 3Ag + 2Ag)^ 
(-C + I6A4 + 9Ag + 4 A 2 ) 

^-B + 4A^ - 3Ag + 2A2^^ 
(C + 16 A 4 - 9Ag + 4 A 2 ) 



= R. 


( 22 ) 


(23) 
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2C - 2Ag - 


(24) 


+ Ag + A2 ®i”(20jjjax) + ® ®max ^max 

4A4 cos(40^^J + 3Ag cos(30^ J + 2A2 cos(20^^^) + B cos 0^^^ = 0 


(25) 

(26) 


-A 4 cos(40^^^) - Ag cos(30^^^) - Ag cos(20^^^) + C cos 0 ^^^ 

+ A4 - Ag + Ag + C = Zj, (27) 

Equations (22) to (27) can be solved using Newton’s method to obtain A 4 , Ag, Ag, 
B, C, and ^n» ^max’ ^b’ h determined from fig- 

ure 21 and then used to calculate the constants in equations (21). .The approximating ana- 
lytic body shape was then plotted and it appears in figure 22. Note that the afterbody in 
this figure is a better approximation to a blimt, spherically capped cone than to a hemi- 
sphere. This condition can be improved by specifying a much larger radius of curvature 
for the base than the true Rg of the actual vehicle, Also, the forebody cone angle and 
the nose radius-of-curvature specification for the analytic body were readjusted to agree 
more closely with the real body. The resulting analytic body shape (fig. 23) has an after- 
body which looks like a hyperbola and which more closely approximates the shape of the 
real configuration as shown in figure 24. 


METRIC COEFFICIENTS 


The metric coefficients of an orthogonal curvilinear coordinate system are of the 
nature of scale factors which give the ratios of differential distances to the differentials 
of the coordinate parameters (ref. 8 ). The equations which are used to define the metric 
coefficients can be derived in the following manner. Consider the orthogonal transforma- 
tion as defined in equations (1) for axisymmetric bodies. Let a 0-curve be defined in space 
as a line of constant r and constant 0 in equations (1). If 7 is the position vector of 
a point P in space, then a tangent vector to the 0-curve is given by 




9Sg d 0 


(28) 


where s^ is the arc length along the 0 -curve. The vector is a unit vector tangent 

6 

to the 0-curve which is redefined as u^ (see fig. 25). Equation (28) can now be rewritten 

U 

as 
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(29) 


® = "e‘‘9 


ds 


n 

where = — - is the length of 9. Therefore the relation for h« can be written as 


de 


h 

^ d0 


dr 


90 



90 90 90 


(30) 


The geometric relationship between 0, and h^ on the 6 -curve is shown in 
figure 25. The final expression for h^, along with h^, and h^ which can be obtained 
in a similar manner, becomes 


^ \ 90 / \ 90 / (90 


N 

(B sinh r - C cosh r) sin ^ ^ sin(n0) 


n=2 


N 

(B cosh r - C sinh r) cos ® ^ cos(n0) 

n=2 


(31) 


^ Ur 


2 \2 r ^ 

(^) (^) ” r + C sinh r) cos ^ ^ cos(n0) 


n=2 

N 


(B sinh r - C cosh r) sin 0 + ^ nAj^e^^ sin(n0) 


n=2 


(32) 


/ v2 /. \2 /ft \2 r ^ 

1 , = ( — 1 + (—1 + | — I = (B cosh r - C sinh r) sin 0 + / 

4> \d4>) \d(f)j \9(^/ r ^2 


sin(n0) 


(33) 


2 2 

We note that hg = h^ in equations (31) and (32). The metric coefficients are all 
positive in a right-handed coordinate system so long as Ug is in the same direction 
as *0, and is in the same direction as R. Therefore, it can be written that 
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TRANSFORMATION OF UNIFORM VELOCITY FIELD 


Note that in figure 20 the direction of increasing r is towards the body and the 
direction of increasing 9 is counterclockwise around the body. The velocity compo- 
nent u is directed along lines of constant r in the direction of increasing 0, The 
velocity component v is directed along lines of constant 9 in the direction of 
increasing r. Let a uniform velocity field of magnitude approach a body in the 

xy-plane at angle of attack a. All angles are measured in standard trigonometric 
fashion from 0 = 0 in a counterclockwise direction. Let be defined as the angle of 
the vector tangent to a line of constant r in the direction of increasing 9. In figure 26 
it is clear that 



From figure 27 and equations (32), (34), and (36), it can be shown that 
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cos p = 


dx 


N 

(B sinh r - C cosh r) sin 0 + ^ sin(n0) 

n=2 


+ dy^ 


( 37 ) 


sin i// = 




dy 

'dx2 + dy2 


N 

(B cosh r - C sinh r) cos ® ^ cos(n0) 

n=2 


TRANSFORMATION OF FLOW EQUATIONS 
AND BOUNDARY CONDITIONS 

An example of the procedure used in applying the generalized 
system can now be presented. Consider the problem of calculating 
rounding the axisymmetric body described in figure 23 as this body 
atmosphere at hypersonic speed. 

The body shape defined in figure 23 was generated with the following parameters 
substituted into equations (22) to (27): 

R„ = 1.10 9^ = 53° 

= 3.55 = Zc(^c»^n>^max) 

Equations (22) to (27) were then solved using a Newton-Raphson iteration scheme 
to obtain 


orthogonal coordinate 
the flow field sur- 
moves through the 


A 4 = 0.068805390 
A 3 = -0.071121784 
Ag = -0.047149036 

B = 2.286877227 
C = 1.703878216 



These constants completely specify the coordinate system which appears in fig- 
ure 23, and are also substituted into equations (31) to (33) to obtain the metric coefficients. 

In order to simplify the presentation here, it is assumed that the body is at 0° angle 
of attack to the incoming flow and that the fluid behaves as a perfect gas everywhere. 
Furthermore, it will be assumed that there is no radiative energy transfer nor any inter- 
nal heat generation in the flow field. 

By starting with the governing equations written for a general curvilinear orthogonal 
coordinate system (refs. 2 and 9), the following equations can be obtained: . 


Continuity 


X • . z 


Xj- momentum 


r~ v”*' 

±(hjh2h3pa.j] + +pu,i2 - 

9 r n / 2 N^^2 

93^ ^1^3(P‘^l"2 ■ ”^2l) . " ■ '^22)^ 

" L. -J 1 

9hg 9h, 

"^2(P - ^33)^ + h(P^1^2 - ^2l)i^ = ° 


X 2 -momentum 






(38) 


(39) 


(40) 
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Energy 


* iT ■ "l"ll - “2^21 + 'll) 


9Xr 


- ^1^12 - ^2^22 ^^2) 


= 0 


Equation of state 


y “ 1 ^ 

P =- P 

y 




where 


Tfi = 2jLi(^ ^ + ^/-i-/h2h„w,) + — (hihgWg 

9Xj hjhg 9x2y • hj^hghg^SxjV 2 3 1 ^ 8 x 2 ^ ^ ^ ^ 


,, = 2m(-1.^ + — /— (h,h,is.) + — (h.h.u. 


Too - 


T 33 = 2n 


Wf ahg ^2 ahgX ^ /g 


hghfaxi h^h^hgl^axi 


— (h2h3Wi) + — (hjh3U.2 


T = T = u 1 ±-l^ + !^1 

12 21 Uiax^Ug) h2 9x2L 


q. = -^J- 
I p h. 9x. 
^r 1 1 


B-|K^*“2^) 


( 41 ) 


(42) 


(43) 


(44) 


(45) 


(46) 


(47) 
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In equations (38) to (47), the Xj coordinate is defined as B and the Xg coordi- 
nate as r. 

The metric coefficients can now be calculated from equations (31) to (33) and sub- 
stituted into the governing differential equations with 


hi = hg = h 

(48) 

h2 = h^ = h 

(49) 

^3 " % 

(50) 

3 _ ^ 

(51) 

dXj^ 90 


9^9 

(52) 

3X2 ■ 8r 


Wj^ = u 

(53) 

Wg = V 

(54) 


(Derivatives with respect to 0 are equal to 0 since this is an axisymmetric flow.) 

In figure 23 the stagnation streamline for 0° angle of attack coincides with the lines 
of constant 9 for 0 = 0 and 0 = tt. For 0 = 0 or 0 = tt, it can be seen from equa- 

3h , 

tions (31) to (33) that h^ = 0, 86 ~ ^ stagnation line. The limiting 

form of the governing equations on the stagnation streamline can now be obtained by 
differentiating equations (38) to (41) with respect to 0. Thus, for 0 = 0 or 6 = n, 


Continuity 






at 




a0 


'(55) 
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0 - momentum 


u = 0 


r-momentum 


at 




uv 



r\ 

r / 0 n1 

^12) 

a 

Jr 

h(p +pv^ - ^ 22 ) 





39 dr 
39 



+ T, 


33 



( 56 ) 


(57) 


Energy 


3_ 

3t 


h^pH 


.2-L 

39 


h/puH - UTj^j - VT2J + 


ar 


hfpvH - VT22 + 


V 


^2) 


a^h 




+ h(pvH-VT22+q2)-^=0 

-± 


39 


where on the stagnation line 


(58) 


^12 = ^21 = ° 


(59) 


. , = 2Ji ^ + v_ X |±(hu) + h^ + ^(hv) + hvii-'" 

11 ^lh30. h2 9rJ h2\39 39 3r 9h 


(60) 
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(61) 


(62) 

The boundary conditions on the body are 

No slip 

u(0) = 0 

v(e) = 0 
Adiabatic wall 


^22 - 9r) h2 


f 8 ^ 

l.(hu) + h^ + i-(hv) + hv-^^^ 
80 80 8r 9h. 


A 

- 9„I 1 9u V 96 9r 

'33 = 2^\hi6 

80 


d\, , 

± 

80 8r 


+ -^1 — (hu) + h— + — (hv) + hv 
h2\80 80 8r 9h 


80 


8T 

8r 


(0) = O 


The boundary conditions at 0 = 0 (the afterbody stagnation line) and 0 = v (the 
forebody stagnation line) are 

u(r) = 0 


^(r) = 0 

80 


^(r) = 0 

80 


^(r) = 0 

80 


M(r) = 0 

80 
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Outflow boundary conditions must be applied at r = r^. The manner in which the 
boundary conditions can be specified or approximated is a question which must be exam- 
ined. However, it should be noted that due to the exponential stretching of the mesh as 
the magnitude of r increases, nominal values of can move the outflow boundary 
far away from the body. For example, if some approximate outflow bovuidary conditions 
can be calculated at a distance of 100 body lengths downstream, then r^ can be calcu- 
lated for this problem using equations (4) as r^ = -5.181. 

It is also possible to map the outflow boundary to infinity and retain orthogonality 
by using an additional coordinate transformation. For example, the coordinate 
transformation 

r = a In 77 


will map 


-oo<r<OtoO<77gl 
where a is an arbitrary constant. 

An explicit time -dependent scheme with appropriate initial conditions can now be 
used to program the governing equations in the rectangular domain of figure 23. (The 
details of any particular solution technique and the computational problems associated 
with it are not the subject of this paper.) The equations can then be integrated forward 
in time until a steady state is reached and a final flow-field solution is obtained. 

CONCLUDING REMARKS 

A transformation procedure has been presented which has the capability of making 
analytic approximations to axisymmetric and two-dimensional body shapes of interest to 
aerodynamicists. The transformation from a rectangular domain in the 0r-plane results 
in an orthogonal coordinate system in the xy -plane which surrounds the specified body 
shape. Configurations which can be modeled by this procedure include spheres, axi- 
symmetric ellipses, spherically capped cones, flat-rfaced cylinders with rounded corners, 
circular disks, and planetary probe -vehicle shapes. Good approximations to spherically 
capped cones are obtained for bodies having a ratio of nose radius to base radius of 1/2 
and cone angles varying from 55° to 60°. 

This transformation procedure should be well adapted to vector-processing com- 
puter systems because the body boundary conditions can be easily calculated along one 
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boundary of a computational plane. Also there is no need to divide the computational 
plane into a forebody flow region and an afterbody flow region. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
February 16, 1977 
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r-rb(0) 

'rs(e)-r|,(e) 



Figure 2.- Polar transformation for blunt -body calculations. 
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(a) Spherically capped cone with hemisphere afterbody. 



Coordinates of A 


Coordinates of B 



(c) Body with concavity. 


Figure 3.- Body-oriented coordinate system. 
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(a) C = 1. 

Figure 6.- Family of ellipses for B = 1. 
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(b) Ag = 0.5. 
Figure 7.- Continued. 
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Figure 8.- Family of triangular shapes for 
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Figure 8.- Continued. 
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Figure 8.- Continued. 
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Figure 8.- Concluded. 
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Figure 9.- Family of rectangular shapes for Ao = 1/9, B = 1, and 
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Figure 10.- Family of rectangular shapes for 
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Figure 10.- Concluded. 
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Figure 11.- Family of pentagonal shapes for N = 4. 
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Figure 12.- Family of hexagonal shapes for N = 5. 
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Figure 13.- Geometry of spherically capped cone. 
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Analytic Body 
Hyperbola 
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Figure 18.- 30° sphere-cone approximation. 
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R„- 1.1875 
Rb - 2.25 

Sr -60° 
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Analytic Approximation 
Planetary Probe Vehicle 



Figure 22.- Planetary probe-vehicle approximation. 
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Analytic Approximation 
Planetary Probe Vehicle 



Figure 24.- Adjusted planetary probe-vehicle approximation. 
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Direction of Increasing 9 
Positive u 



Figure 26.- Free -stream velocity components along line of constant 
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